Small strange stars and marginally stable orbit in Newtonian theory 



Julian Leszek Zdunik 

N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warszawa, Poland, 

e-mail: jlz@camk.edu.pl 

Eric Gourgoulhon 

Departement d'Astrophysique Relativiste et de Cosmologie - UMR 8629 du CNRS, Observatoire de Paris, 
F-92195 Meudon Cedex, France, e-mail: Eric.Gourgoulhon@obspm.fr 
(15 December 2000) 



It is shown that for very rapidly rotating low mass strange 
stars the marginally stable orbit is located above the stellar 
surface. This effect is explained by the very important role of 
the oblateness of the rotating strange star. The comparison 
with some "academic" examples is presented. This feature is 
purely Newtonian in its nature and has nothing to do with 
relativistic marginally stable orbit. The effect is very large 
and cannot be treated in a perturbative way. It seems that 
strange stars as a very dense self-bound objects are the only 
possibility in Nature to represent these toy models. 
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Recently marginally stable orbits have been exten- 
sively studied as a crucial input of the most popu- 
lar model explaining kilohertz quasi-periodic oscillations 
(QPOs) in X-ray binaries (see [l| for a review). Besides 
neutron stars as a possible source of these QPOs hy- 
pothetical strange stars have been considered. Strange 
stars, objects built of the mixture of u, d, s quarks, 
in contrast to neutron stars have no minimum mass; 
their size extends continuously from that of large nuclei 
(strangelcts) up to the maximum mass configurations of 
about 2 solar masses (for a recent review, see ||]). It has 
been showed ||^] that the marginally stable orbit is lo- 
cated above the surface of the strange star for relatively 
high mass even at very high rotation rates, which is the 
main difference with respect to neutron stars |^|. It was 
suggested pj that, for configurations close to the Keple- 
rian one, defined as a mass-shedding limit at the equator, 
the very large deformation of the strange star leads to the 
significant increase of the radius of the marginally stable 
orbit, resulting in a gap between this orbit and the stellar 
surface. Of course in these considerations, performed for 
relatively massive (M > M Q ) stars, the GR effects seem 
to play a dominant role. General relativity predicts that 
the radius of the marginally stable orbit for nonrotating 
star of mass M is equal to « 8.8M/M© km, which is com- 
parable to the radius of a strange star. It is well known 
that quadrupole moment of the star could also destabi- 
lize the orbital motion [||, but it seemed that this effect 
is not very large in real physical situations, at least as far 
as neutron stars are concerned. We show that for objects 
such as self-bound strange stars the deformations due to 
the rotation can be so large that the change in Newto- 



nian gravitational potential could result in the existence 
of a marginally stable orbit, even when GR corrections 
are completely negligible. 

The existence of a marginally stable orbit corresponds 
to the conditions dV/dr — and d 2 V/dr 2 — where 
V(r) is the effective potential (for the radial part of mo- 
tion) for a particle orbiting the central body at a fixed 
angular momentum per unit mass I. The resulting equa- 
tions are discussed in standard textbooks in the case of 
the Schwarzschild geometry (see e.g. In first approx- 
imation, the source of the marginally stable orbit is the 
term in V proportional to 1/r 3 with the same sign as 
the gravitational potential (i.e. opposite sign to centrifu- 
gal force). In general relativity it is a natural result of 
constants of geodesic motion on circular orbits. 

In Newtonian theory this term is generated by the 
quadrupole moment of the star Q\ in the equatorial plane 
the gravitational force per unit mass of the orbiting par- 
ticle is 
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With the above convention Q > for an oblate star and 
Q < for a prolate one. The relativistic definition of Q 
can be found in Ref. |@|. 

One can take into account the effect of the oblateness 
and of the angular momentum of the star using standard 
methods, i.e. via expansions of the metric functions in 
powers of 1/r 2 . In the framework of GR this leads to the 
following formula for marginally stable orbit: 



la+yr+gGR/e)- 



3^ g - ■ _ : ( 2 ) 

where q GR = Qc i /G 2 M 3 , j = Jc/GM 2 , R s = 2GM/c 2 , 
J being the total angular momentum and we linearized 
the dependence in j. In the case of Q = Eq. (||) reduces 
to the formula derived by Kluzniak and Wagoner ||. 
In nonrelativistic limit we have 



Ru 



R e 



Rl 



2Q\ 
3MJ 



3Q 
2M 



(3) 



In Newtonian theory the marginally stable orbit can be 
determined exactly for some toy models. We will present 
here results for one- and two-dimensional cases. 
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In one dimension we consider a bar with given constant 
linear density pi, i.e. mass of the bar is M = 2piR, 
where R is the "radius" of the bar equal to the half of 
its length. The gravitational attraction on the line which 
coincides with the bar at the distance r from the bar 
center is equal to f g = —GM/(r 2 — R 2 ). Of course we 
have to assume that this field does not change during 
orbital motion of the particle, i.e. the bar rotates with 
the same angular velocity as the orbiting particle. The 
radius of the marginally stable orbit is then 
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R^f = V3R -1.73R. 
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FIG. 1. Effective potential of the particle on the circular 
orbit around the disk when the square of angular momen- 
tum I 2 is equal to 1.6, 1.8 and 2GMR. The solid line corre- 
sponds to the case of marginally stable orbit (the exact value 
is l'i s = 1.784 GMR). 

In two dimensions the corresponding body is an in- 
finitesimally thin disk of constant surface density p s with 
mass M = Trp s R 2 . The gravitational force in the plane 
of the disk at the distance r is then given by the hyper- 
geometric function 2 F\ 
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where q = 2rR/(r 2 + R 2 ). The form of the effective 
potential [the integral of Eq. (|5|) plus a centrifugal term] 
results in the marginally stable orbit located at the radius 
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The angular momentum of the particle at this orbit is 
then equal to I ~ VX784GM R (per particle mass). 

It should be mentioned that in these examples expan- 
sions of the gravitational potential in powers of (R/r) 2 
would be completely unjustified. For example in the case 
of a disk the expansion of Eq. (||) leads to the formula: 
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Using only the second term (quadrupole order) we obtain 
the radius of marginally stable orbit 



i.e more than twice smaller than the exact value. 
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FIG. 2. Equatorial radius of the strange star R (solid), 
radius of the marginally stable orbit _R ms (dotted), quadrupole 
approximation of J? ms - i?m" ad (dashed), and radius of the 
marginally stable orbit if the rotating object would be a disk 
of the same radius i?^ 11 . The star has a mass 0.005Mq. The 
value of R^t at J = corresponds to the relativistic radius 
of the marginally stable orbit for nonrotating star which in 
this case is equal to 0.044 km. The relativistic term due to 
the dependence on j (Eq. ^) is negligible. 

In previous papers ||f| the marginally stable orbit for 
rotating strange stars have been determined in the frame- 
work of General Relativity and the authors concluded 
that the main difference with respect to neutron stars is 
the gap between this orbit and the stellar surface for high 
rotation rates. 

Rotating strange stars are very oblate for large rota- 
tion rates and the Keplerian configurations for a sequence 
of stars with given baryon number is achieved not due to 
the increase of the rotation frequency v T0 t, but due to 
the deformation of the star Q. It was suggested that 
this large oblateness is the true reason which cause the 
radius of the marginally stable orbit to be above the stel- 
lar surface. 

In this paper we try to prove this statement by consid- 
ering strange stars for which relativistic effects arc very 
small. Thus we will focus on the results for small strange 
stars where the oblateness is crucial for the properties 
of these objects and in particular, is the source of the 
marginally stable orbit. The effects of General Relativ- 
ity, which are very important in massive dense stars, are 
here negligible. However presented results has been ob- 
tained in the framework of GR using the numerical code 
developed by Gourgoulhon et al. p0| . 

We calculate rotating bare strange stars (without 
crust) using the strange matter model considered as a 
basic one in previous papers (and called SQM1 in Q). 

In Fig. we present the equatorial radius of the star 
R cq with fixed baryon mass M — 0.005 M Q and the ra- 
dius of marginally stable orbit. We choose such a small 
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object in order to extract Newtonian effects - which for 
massive stars are contaminated (and for slowly rotating 
stars dominated) by relativistic ones. We see that for 
high rotation rates there exist a gap between stellar sur- 
face and the marginally stable orbit. The curve Rf^l k is 
defined by Eq. (g) with R = R cq . The curve i?Sf d corre- 
sponds to the quadrupole approximation of i? ms given by 
3) with Q calculated in the framework of GR [Eq. (7) 
We see that this quadrupole approximation is 
wrong and as in the case of disk the resulting radius is 
about twice smaller than the true value. This conclu- 
sion confirms the results of Shibata and Sasaki ||, who 
found necessary to include higher multipole moments in 
the marginally stable orbit determination in the case of 
polytropic EOS. They concluded that the quadrupole 
moment is as important as the spin. In our model the 
effect of the rotational deformation of the star is larger 
because the matter is self-bound at very high density and 
the stars are more oblate. Due to these two reasons for 
fast rotating strange stars we are "closer" to the disk ap- 
proximation than to the expansion up to the quadrupole 
term (see Fig. ^). The value of the quadrupole moment 
for small strange stars close to the Keplerian configu- 
ration is Q m ax = 0.1 7 MR 2 which is quite close to the 
"disk" value - 0.25MR 2 . 



The behavior of R cq and R ms for a star rotating 
faster and faster, presented in Fig. |5[ is typical for 
small mass stars self-bound at high density. For higher 
masses (larger than 1.12 M© in this model) the radius 
of marginally stable orbit is larger than R at 2/ ro t = 
and for low rotation rates decreases due to the relativis- 
tic effects and then increases due to the oblateness (see 
H for figures and discussion). Of course the point at 
which -Rms^rot) gets larger than the equatorial radius 
R eq depend on the mass of the star under consideration. 

In the next figures we present the regions of stellar 
parameters and rotation rates for which marginally stable 
orbit is located above the stellar surface. 

In Fig. m we show the boundaries resulting from the 



Keplerian frequency and from the condition R c 
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We see that in M — R plane the region with R ms < R eq 
looks rather small and even for low-mass stars there ex- 
ist relatively large area in which R cq < R ms . It is also 
clearly visible that the stars with masses between 1.12 
and 1.28 M Q enter the region with R ms < R cq during 
slowing down or acceleration although for high and low 
rotation rates R ms > R cq . These stars for fixed baryon 
mass are represented by nearly horizontal (vertical) lines 
in top (bottom) panel of Fig. ^ ( "nearly" because gravi- 
tational mass increases a little due to the rotation). 
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FIG. 3. The regions where equatorial radius of the 
marginally stable orbit is smaller than stellar radius 
(hatched). Solid curve is defined by the equality R eq = -Rms. 
The condition R cq < R ms is fulfilled between solid and 
dotted curve (corresponding to the Keplerian frequency). 
The dashed curve represents nonrotating configurations (top 
panel). In hatched area the ISCO is defined by the stellar 
surface. The region for 1.12 < M/M Q < 1.28 corresponds 
to the configurations for which for given mass _R ms decreases 
with Vrot due to dependence on j, then crosses the radius and 
for higher u Iot reappears. 
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FIG. 4. The region in the i/ ms vs M plane where the ra- 
dius of the marginally stable orbit is smaller than the stellar 
radius (hatched area). In the double hatched region there ex- 
ist configurations with R cq < R ma (slowly rotating) as well 



(see text). The solid curve is defined by the 
= -Rms, the dashed line is the frequency at the 
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ISCO for nonrotating star, which for M < 1.12Mq coincides 
with stellar radius and for higher masses is proportional to 
1/M. The long dashed curve is the maximum i/ ms for higher 
mass strange stars when R eq < R ms everywhere. The dot- 
ted curve is the minimum value of i/ ms achieved for Keplerian 
configurations. 

In Fig |] we present limits on the frequency of 
marginally stable orbit taking into account the condi- 
tions i? e q < -Rms which in prin ciple can be supported by 
the observations of QPOs |pT[. 
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There exists a quite large region in the v ms -M plane, 
where this condition is fulfilled. It should be mentioned 
that because z^ ms in non-monotonic function of J for a 
star with given baryon mass, one region in v ms -M plane 
is not uniquely defined. Namely in double hatched area 
there exist stellar configurations with R cq < R ms (slowly 
rotating) as well as stars for which R cq > R ms . For non- 
rotating stars we have here R eq < i? ms but if such a star 
rotates faster and faster it moves upwards, turns back at 
the solid line and moves downwards having R cq > R ms 
as long as it do not leave the region between solid lines. 
One has similar situation (non-uniqueness) in the case of 
high mass stars (M > 1.28M Q ), which turn back at long- 
dashed line but the condition R cq < i? ms is fulfilled all 
the time. In non hatched region all stars have R eq < Rms- 
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FIG. 5. The ratio T/W (solid) and the quadmpole moment 
Q (dashed) as a function of the oblateness of the star along 
the sequence of rotating strange stars for which R aq — i? ms 
(solid curves in Figs Point (0,0) describes nonrotating 
configuration. The last points (maxima) correspond to the 
smallest stars. The Keplerian limits are located far above and 
to the right of this graph (T/W ~ 0.25, oblateness ~ 0.63) 

From Fig. || we see that the oblateness of the star, 
the quadrupole moment Q and the ratio of the kinetic to 
the potential energy T/W for the rotating stellar config- 
urations defined by the equality R eq = R ms are highly 
correlated and almost linearly dependent. Although the 
frequency at which this condition is fulfilled is approx- 
imately constant for M < 1 M Q (0.9-1 kHz, Fig |) the 
stars with lower mass are much more oblate at this rota- 
tion rate, having proportionally larger T/W and Q. For 
these configurations large value of R ms results entirely 
from oblateness of the star. Albeit the role of the oblate- 
ness for the location of the marginally stable orbit is most 
pronouncing in the case of small strange stars close to the 
Keplerian velocity, also for massive stars rotating with 
relatively high frequency the deformation of the star is 
crucial for the existence of the marginally stable orbit. 

In this paper we showed that the gap between inner- 
most stable orbit and the surface of the star can be pro- 
duced by the oblateness of the star in the framework 
of Newtonian theory. Due to the compactness and very 



high surface density strange stars are natural candidates 
for this effect to be a real possibility in Nature. This 
is due to two reasons. First: the stars are very oblate, 
and the gravitational force is far from the attraction by 
a point mass. Second: the stars are self-bound at very 
high density and in practice an oblate small strange star 
is quite close to a constant density disk. In contrast to 
bare strange stars, for oblate neutron stars (and also for 
strange stars with crust) the low-density region is rela- 
tively large making the radius of the star larger than R ms 
(or in other words the radius of corresponding "disk" is 
much smaller than the stellar radius). 
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